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1.  Introduction  and  Summary 

Readiness  evaluation  is  one  of  the  most  important  problem  areas  in 
the  study  of  complex  military  systems.  Such  studies  usually  encompass  a 
large  number  of  measurements  on  the  performance  of  the  various  subsystems  and 
then  attempt  to  construct  reasonable  models  that  relate  the  evaluation  indices 
of  the  subsystems  in  a  meaningful  functional  manner.  This  is  indeed  very 
often  a  formidable  task  due  to  the  complexity  and  multiplicity  of  variables 
and  functions.  However,  it  is  often  the  case  that  many  of  the  measured  vari¬ 
ables  correlate  with  each  other.  These  intercorrelations  reveal  that  variables 
contain  some  information  on  each  other.  Accordingly,  if  these  intercorrelations 
can  be  utilized  in  a  manner  that  allows  considerable  reduction  in  the  number 
of  factors  to  be  considered,  without  much  loss  in  the  information  in  the  orig¬ 
inal  data,  a  significant  step  can  be  taken  towards  simplification  of  the  prob¬ 
lem.  The  present  paper  applies  several  well-known  multivariate  statistical 
methods  to  attain  this  goal.  The  main  objective  of  the  present  paper  is  to 
discuss  what  some  of  the  available  multivariate  statistical  methods  can  attain 
and  to  show  that  such  methods  can  be  easily  implemented  by  utilizing  appropriate 
computer  packages.  In  particular,  we  discuss  the  methods  of  principal  and 
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rotated  factor  analysis  and  we  apply  these  methods  to  simulated  data  on  21 
operational  readiness  variables  related  to  Navy  destroyers.  The  variables 
and  the  corresponding  parameters  were  taken  from  the  Institute  of  Naval  Studies 
study  [9].  This  study  analyzes  actual  data  collected  over  several  years  on  83 
destroyers.  It  extends  to  various  aspects  of  the  readiness  problem  and  relates 
operational  readiness  to  material  readiness.  As  mentioned  in  our  recent  survey 
paper  [5J,  we  believe  that  study  (9]  is  of  fundamental  importance.  It  employs 
a  variety  of  multivariate  statistical  procedures  in  a  penetrating  manner  and 
provides  a  sound  analysis  of  complicated  problems.  Our  intention  is  not  to 
duplicate  that  study  but  to  provide  an  exposition  on  the  application  of  the 
multivariate  methods  mentioned  above.  We  have  chosen  to  create  data  sets  by 
simulation  and  not  to  use  actual  data  since  in  this  way  we  can  generate  data 
following  multivariate  normal  distributions  having  specific  structures.  Thus, 
by  applying  the  multivariate  methods  on  different  sets  of  simulated  data  we  can 
illustrate  the  strength  of  the  methods  and  what  can  actually  be  achieved.  We 
will  show  that  the  systems  (destroyers)  in  this  example  can  be  classified  accord¬ 
ing  to  the  values  of  two  or  three  factor  scores,  which  relate  all  21  variables 
in  an  orthogonal  fashion.  The  factor  scores  can  be  graphed  and  their  periodic 
determination  can  provide  important  follow-up  on  the  state  of  readiness.  Sta¬ 
tistical  control  charts  can  be  devised  to  provide  early  detection  of  deteriora¬ 
tion  in  the  state  of  readiness.  Similarly,  if  the  data  consists  of  a  mixture 
of  two  or  more  samples  from  different  multivariate  populations,  the  plotting  of 
factor  scores  obtained  by  a  factor  analysis  of  the  whole  data  set  can  reveal  the 
existence  of  different  clusters.  These  ideas  will  be  demonstrated  in  the  pre¬ 
sent  paper.  We  start  in  Section  2  with  a  description  of  the  simulations  and  the 
structure  of  the  data  sets.  Section  3  is  devoted  to  principal  and  rotated  factor 
analysis.  In  Section  4  we  discuss  the  application  of  factor  analysis  to  detect¬ 
ing  changes  in  the  state  of  readiness  of  systems.  The  mathematical  development 
is  presented  in  appendixes. 

The  implementation  of  multivariate  analysis  of  the  type  discussed  in  the 
present  paper  previously  required  development  of  computer  programs  and  systems 
for  data  storage  and  analysis  and  this  hindered  its  growth  and  acceptance.  Pre¬ 
sently  there  are  several  statistical  computer  packages  which  are  available  in 
many  computing  centers  and  which  are  easily  implement able.  The  well-known 
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Biomedical  Computer  Programs  (P-Series,  1977)  [7],  Statistical  Packages  for 
the  Social  Sciences  (SPSS)  [11],  and  the  Statistical  Analysis  System  (SAS) 
[3J  are  the  most  suitable  for  our  purposes.  Fortran  programs  for  factor 
analysis  and  other  related  multivariate  techniques  can  be  found  in  Cooley 
and  Lohnes  [6],  Overall  and  Klett  [10]  and  more  specifically  for  the  methods 
used  in  the  present  paper  in  our  report  [4].  In  the  present  paper  we  apply 
the  SAS  procedures.  In  the  appendixes  we  discuss  and  also  present  the  SAS 
programs  that  we  have  used. 

2.  Simulating  Data  Sets 

In  the  present  study  we  construct  data  sets  on  the  basis  of  the 
operational  readiness  indices,  ORI,  of  the  following  21  variables. 


V1 

Ship  control 

SKC 

V2 

Navigation 

NAV 

V3 

Surface  operations  -  CIC 
(Combat  Information  Center) 

SOPS 

v4 

Battle  communications 

BATC 

V5 

Surface  gunnery  (non-firing) 

SGUN 

V6 

AAW  (Anti-air  Warfare)  -  CIC 

AAWC 

V7 

AAW  -  Weapons  Control 

AAWN 

V8 

Engineering 

ENG 

V9 

Setting  material  condition 

SMC 

V10 

Damage  control 

DC 

V11 

NBC  (Nuclear,  biological,  and  chemical) 

NBC 

V12 

Low-visibility  piloting 

LVP 

V13 

CIC  -  Assistance  in  piloting 

CICAP 

V14 

CIC  -  Assistance  in  ASW 
(Anti-submarine  warfare) 

CICASW 
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V15 

ECM  (Electronic  countermeasures) 

ECM 

V16 

Modified  full -power  run 

BFPR 

V17 

Surface  firing 

SFIR 

V18 

AA  firing 

AAF 

v19 

Gunfire  support 

GUNS 

V20 

Communications 

COMM 

V21 

ASW  operations 

ASW 

The  raw  scores  obtained  on  these  variables  by  the  83  ships  during 
training  can  be  obtained  in  the  Institute  of  Naval  Studies  [91.  We  consider 
rather  the  ORl's  which  are  indices  obtained  from  the  raw  scores  by  the 
transformation 

ORI  =5+2  (NSCORE)  (2.1) 


where  NSCORE  denotes  the  standard  normal  fractile  corresponding  to  the 

percentile  point  of  the  raw  score.  More  precisely,  if  x...  <  ...  <  x. 

(1)  =  =  (n) 

is  the  order  statistic  of  a  sample  of  n  observations  on  a  variable  x  , 


the  NSCORE  corresponding  to  x 


(i) 


is 


'(i) 


* 1  U4  ’ 


t  •  i,. 


where  I’(z)  is  the  standard  normal  C.D.F.  Theoretically  the  ORI  values,  of 
each  of  the  variables  v^  (i*l,...,21)  ,  are  normally  distributed  with  mean 

El v ^ }  =  5  (i  =  l,...2l)  and  variance  Variv^}  =  4  (i=l, — 21).  In  addition, 

the  ORI  variables,  v^  ,  are  not  independent.  We  assume  that  the  vector 
T 

X  =  ^VL’'**’V21^  lias  3  mu^t^normal  distribution  N(j^,^)  ,  with  mean  vector 

T 

jji  =  5^  ,  where  ^  =  [l,...,l]  ,  and  covariance  matrix  ^  =  4^  where 

denotes  the  matrix  of  intercorrelations  among  the  21  variables  v^^  .  For  the 
purpose  of  simulating  data  sets  we  have  used  the  matrix  ^  given  in  the 
Institute  of  Naval  Studies  [9]  and  presented  here  in  Table  1.  The  simulation 
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TABLE  l 

I NTE  RCORRE I  .AT I ONS  BETWEEN 
ORI  VARIABLES  v,  t«  v„, 


t  ;  1. 00 
3.  06 

2  0.17 
0.26 

3  0.27 
0.26 

4  0.  04 
0.  03 

5  0.33 
0.  04 

6  0.31 
0.  15 

7  0.  45 
0.  11 

8  0.  10 

0.  01 

V  0.  06 
0.  04 

10  0.03 
0.  07 

11  0.  03 
0.  11 

12  0.32 
0.  06 

13  0.36 
0.  12 

14  0 .  06 
1. 00 

15  0.17 
0.  16 

16  0. 16 
0.23 

17  0. 07 

0.  04 

1«  0. 12 
|  3.03 

19  !  0.03 

0.  06 

20  0.  01 

0.  06 

21  3.10 

I  0.  05 


0.  17 
0.  17 
1.00 
0.  01 
0.22 
0.  01 
0.  06 
0.  13 
0.  13 
0.  04 
0.  03 
0.  06 
0.  14 
0.  17 
0.22 
0.  03 
0.  01 
0.  14 
0.20 
0.20 
0.  04 
0.  03 
0.  12 
0.  13 
0.  14 
0.  04 
0. 26 
0.  16 
0.  01 
1.  00 
0.  07 
0.  13 
0.  02 
0.  1  4 
0.  03 
0.  14 
0.  03 
0.  02 
0.  03 
0.  34 
0.  13 
3.31 


0.27 
0.  16 
0.22 
0.  07 
1. 00 
0.  07 
0 .  06 
0.25 
0.  13 
0.25 
0.  03 
0.  17 
0.  14 
0.  14 
0.22 
0.  10 
0.  01 
0.21 
0.20 
0.  17 
0.  04 
0.  03 
0.  12 
0.25 
0.  14 
0.20 
0.26 
0.23 
0.  01 
0.  13 
0.  07 
1.00 
0.  02 
0.  01 
0.  03 
0.  14 
0.  03 
0.  03 
0.  03 
0.23 
0. 13 
3.  13 


0.  04 
0.  07 
0.  06 
0.  02 
0.  06 
0.  02 
1.00 
0.  03 
0.23 
0.  1  0 
0.23 
0.  12 
0.  13 
0.  01 
0.  1 1 
0.  17 
0.  13 
0.  13 
0.  1  0 
0.  03 
0.  03 
0.  07 
0.  12 
0.  07 
0.  35 
0.  1 1 
0.  03 
0.  04 
0.  13 
0.  14 
0.25 
0.  01 
0.  03 
1 .  00 
0.  10 
0.  00 
0.  03 
0.  05 
0.  04 
0.24 
0.  06 
0.2  3 


0. 33 

0.  12 

0.  13 

0 .  03 
0.  13 
0.  03 
0.23 
0.10 
1 . 00 
0.20 
0.34 
0. 13 
0.  05 
0.31 
0.  03 
0.23 
0.  04 
0.  05 
0.  01 
0.  02 
0.  12 
0.  04 
0.  07 
0.  07 
0.  03 
0.  03 
0.  04 
0.  03 
0.  04 
0.  14 
0.25 
0.  14 
0.  10 
0.  00 
0.20 
1.00 
0.  03 
0.11 
0.  13 
0.  04 
0.  02 
0.  12 


0.  31 
0.  03 
0.  03 
0.  03 
0.  03 
0 .  03 
0.23 
0 .  03 
0.34 
0.  03 
1 . 00 
0.  04 
0.23 
0.  01 
0.  13 
0.  12 
0.  13 
0.  07 
0.  03 
0.  02 
0.  05 
0.  13 
0.  03 
0.  01 
0.22 
0.20 
0.  15 
0.  06 
0.  06 
0.  02 
0.  17 
0.  03 
0.  12 
0.  05 
0.  13 
0.11 
0.  04 
1. 00 
0.  03 
0.  12 
0.  10 
0.  03 


0.  45 
0.  01 
0.  14 
0.  03 
0.14 
0.  03 
0.  13 
0.  04 
0.  05 
0.  13 
0.  23 
0.  03 
1 . 00 
0.  12 
0.  17 
0.  07 
0.  17 
0.  13 
0.  03 
0.23 
0.  03 
0.  06 
0.  14 
0.23 
0.  17 
0.  17 
0.  11 
0.  06 
0.  17 
0.34 
0.  14 
0.23 
0.  01 
0.24 
0. 31 
0.  04 
0.  01 
0.  12 
0.  12 
1.  00 
0.  03 
0.  03 


0.  10 
0.  1  0 
0.22 
0.  13 
0.22 
0.13 
0.  11 
0.  06 
0.  03 
0.  02 
0.  13 
0.  I  0 
0.  17 
0.  03 
1 .  00 
0.  03 
0.  03 
0.  13 
0.10 
0.  13 
0.  00 
0.  13 
0. 21 
0.21 
0.  40 
0.  05 
0.  01 
0.  05 
0.  03 
0.51 
0.  10 
0.  13 
0.  17 
0.23 
0.23 
0.  12 
0.  12 
0.  03 
0.  07 
0.  03 
0.  03 
1.  00 


0.  06  0.  03  0.  03  0.  32  0.  36 


0.  04  0.  12  0.  14 


0.20  0.  04 


0.  14 


0 .13 


0.  04 


0.  03  0.  12  0.  35  i 


0.  12  0.  07  0.  03 


0.  13  0.  03  0.  05  0.  03  0.22 


0.  03  0.  03  0.  14 


0.  03  0.  10  0.  00 


0.40 


1.  00  0.  13  0.  13  0.  07  0.  14 


0.  13  1.00  0.63  0.  06  0.04 


0.  13  0.63  1.00 


0.  06 


0.  02 


1.00  0.  00 


0.  14  0.  04  0.  02  0.  00  1 .  00 


0.  04 


0.  14  0.20 


0.13  0.03 


0.  05 


0.  06  0.  12 


0.13  0.04 


0.  03  0.25  0.20 


0 .  04 


0.  02  0.  13 


0.13  0.23 


0.  13 


0.23 


0.  13 


O.'tlS 


was  performed  according  to  an  algorithm  described  in  Appendix  I.  It  is  based 
on  simulating  independent  standard  normal  variates,  z  ,  and  transforming  them 
to  corresponding  variates  (i-l,...,21)  by  a  transformation  involving 

th*'  eigenvalues  and  eigenvectors  of  the  correlation  matrix.  An  SAS  program 
for  such  simulation  is  given  in  Appendix  III.  In  Table  2  we  present  a  sample 
of  n  =  50  vectors  of  six  variables  (v^,  v,.,  v^,  v^,  v^0,  v.^)  simulated 

according  to  this  program.  The  sample  means,  standard  deviations  (STD  DEV) 
and  intercorrelations  are  provided  in  Table  3.  As  illustrated,  the  sample 
statistics  are  generally  deviating  to  some  extent  (according  to  their  sampling 
distributions)  from  the  parameters  used.  However,  in  actual  cases  the  popula¬ 
tion  parameters  are  unknown  and  the  analysis  must  be  based  exclusively  on  the 
sample  values,  with  the  possible  incorporation  of  some  prior  information,  and 
this  is  what  we  are  doing  here. 

It  should  also  be  remarked  that  the  simulation  is  based  on  the  matrix 
of  intercorrelations  of  the  above  six  variables  only.  This  matrix  is,  how¬ 
ever,  a  submatrix  of  that  given  in  Table  1  and  can  be  obtained  by  reading 
the  appropriate  rows  and  columns.  In  the  course  of  the  present  study  several 
different  data  sets  were  simulated,  employing  the  same  algorithm  with  only 
slight  modifications  from  case  to  case,  as  will  be  explained  later. 


3.  Principal  and  Rotated  Factor  Analysis 

It  is  generally  difficult  to  make  comprehensive  inference  of  multi¬ 
variate  data  without  further  analysis,  due  to  the  large  number  of  inter- 
correlated  variables.  Even  in  the  case  of  only  six  variables  it  will  be  diffi¬ 
cult  to  discriminate  between  "good"  and  "bad"  systems,  just  by  inspecting  the 
data  sets,  or  by  performing  a  univariate  analysis  on  each  variable  separately. 
The  methods  of  multivariate  analysis  are  designed  to  provide  the  needed 
information  in  cases  of  many  correlated  variables.  In  the  present  section 
we  discuss  the  methods  of  principal  and  rotated  factor  analysis,  and  show  how 
they  can  be  applied  to  the  evaluation  of  the  readiness  of  systems.  An  outline 
of  the  theory  is  given  in  Appendix  II.  We  refer  the  reader  for  an  extensive 
development  of  the  theory  and  computer  programs  to  the  books  of  Overall  and 
Klett  [10],  Cooley  and  Lohnes  [5],  Tatsuoka  [12],  and  Van  de  Geer  [13]. 
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50  SIMULATED  VECTORS  OF  SIX  VARIABLES 
ACCORDING  TO  THE  MULTINOMIAL  DISTRIBUTION  N(51,4jy 
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TABLE  3 

SAMPLE  STATISTICS  OF  THE 
SIX  VARIABLES  IN  TABLE  2 


V1 

V5 

V6 

V7 

V12 

V14 

MEANS 

4.860 

4.577 

4.644 

4.728 

4.993 

4.914 

STD  DEV 

1.998 

1.897 

2.180 

2.244 

1.916 

2.468 

CORRELATION  MATRIX 


V1 

V5 

V6 

V7 

V12 

V14 

V1 

1.000000 

0.535377 

0.183915 

0.558162 

0.399037 

0.519481 

V5 

0.535377 

1.000000 

0.158194 

0.362251 

0.238443 

0.128076 

V6 

0.183915 

0.158194 

1.000000 

0.237437 

0.192687 

-.023431 

v7 

0.558162 

0.362251 

0.237437 

1.000000 

0.492125 

0.617677 

V12 

0.399037 

0.238443 

0.192687 

0.492125 

1.000000 

0.239948 

V14 

0.519481 

0.123076 

-.023431 

0.617677 

0.239948 

1.000000 

3.1  Principal  Factor  Analysis 

The  main  objective  of  principal  factor  analysis  is  to  provide  a  small 
number,  m  ,  of  linear  combinations  of  the  original  variables 
v^,...,Vp  (2<m<p)  so  that  (i)  a  large  proportion  of  the  total  variance  of 

the  original  variables  should  be  accounted  for  by  the  m  transformed  vari¬ 
ables,  and  (ii)  the  transformed  variables  should  be  uncorrelated.  It  is 
shown  in  Appendix  II  that  the  solution  of  this  problem  is  obtained  by  deter¬ 
mining  first  the  m  largest  eigenvalues  of  ^  and  the  corresponding  eigen¬ 
vectors;  followed  by  determination  of  factor  scores  for  each  system.  Let 
A^  >  . . .  >  Ap  >  0  be  the  eigenvalues  of  the  p  x  p  correlation  matrix  . 

Since  £  is  positive  definite,  these  eigenvalues  are  all  real  and  positive 
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(with  probability  one).  Moreover  A  j+. . .+Ap  =  p  .  Hence,  choose  m  so 

that  (A^. .  ,+Am)/p  is  "close  enough"  to  1.  This  ratio  is  the  proportion 

of  the  sum  of  variances  of  v}  (i=l, . . . ,p)  that  is  accounted  for  (explained) 

by  the  m  factors.  These  factors  are  constructed  in  the  following  manner. 
Let  Jj>^  (j*l,..,m)  be  the  orthonormal  eigenvector  of  ^  corresponding 

to  (j“l»« . . , m )  .  The  m  factor-score  variables  corresponding  to 
T 

X  =  [v1>...,vp]  are  given  by 


j  =  1 . m  (3.1) 


T 

where  =  [u^,...,up]  is  a  vector  of  standard  scores  corresponding  to 

A 

X  »  u^  =  (v^-v^)/o^  ,  i  =  l,...,p  ,  v^  denotes  the  sample  means  of 

/A 

the  ith  variable  and  designates  its  sample  standard  deviation.  How  much 

statistical  information  available  in  the  original  vectors  of  p  variables  is 
retained  in  the  m  factor  scores  of  the  individuals  in  the  sampl>  ?  To  ans¬ 
wer  this  question  we  introduce  additional  concepts  from  the  theory  of  factor 
analysis . 

Consider  the  matrix  ,  of  order  p  x  m  ,  whose  m  column  vectors 
are  related  to  the  m  largest  eigenvalues  and  their  corresponding  eigen¬ 
vectors,  according  to  the  formula: 


*u> '  xj1/2  *0) 


j  "  1 » •  •  •  jin  •  (3*2) 


This  matrix  is  called  the  factor  pattern  (structure)  matrix.  Obviously, 

I  I  I  1 2  *  Vj  ,  j  =  1,. .. ,m  ;  and  m  =  p  then  ^  ,  or 

P  t 

•  %  ~  ^  ^(j)^  (j)  ’  Thls  is  the  sPectral  decomposition  of  the  correlation 

A  ^  m  ^  -A 

m<p  we  ft. '  £  (j)  “nd  ft.  •  ft  '  ft.  • 


.  matrix  .  If 
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It  is  desirable  to  choose  m  so  that  is  negligible  (or  statistically 

Insignificant).  Tests  of  significance  of  are  available  (see  Cooley  and 

Lohncs  [6,103].  According  to  the  spectral  decomposition  of  ^  ,  the  proportion 
of  the  ith  diagonal  element  of  given  by  the  corresponding  element  of 

is  called  the  communality  of  the  ith  variable.  It  is  determined  by  the 


formula 


m 

E 

J-l 


i  =  l,...,p  .  (3.3) 


One  can  say  that  h^  is  the  proportion  of  the  variance  of  the  ith  variable 
v^  explained  by  the  ra  factors.  In  Table  4  we  present  the  eigenvalues  and 
the  corresponding  eigenvectors  of  the  correlation  matrix  of  Table  3.  These 
eigenvalues  and  eigenvectors  were  obtained  by  employing  a  computer  library 
routine.  On  the  basis  of  these  values,  the  first  three  factor  scores  were 
determined  for  the  50  simulated  vectors  of  Table  2,  according  to  formula  (3.1). 
These  factor  scores  are  presented  in  Table  5. 


TABLE  4 

EIGENVALUES  AND  EIGENVECTORS  OF  THE 
CORRELATION  MATRIX  IN  TABLE  3 


Eigenvalues 

X.  X. 

X„ 

X, 

X, 

X 

1 

2 

3 

4 

5 

6 

2.749232 

1.068542 

0.850006 

0.700426 

0.384321 

0.239473 

Eigenvectors 

^(1)  ^(2) 

^(3) 

fe(4) 

^(5) 

^(6) 

0.503687 

-.035982 

0.266720 

-.115832 

0.707698 

0.399522 

0.360454 

0.292523 

0.749462 

0.026661 

-.329692 

-.336753 

0.192362 

0.755713 

-.410358 

-.440037 

0.072605 

-.156839 

0.514023 

-.118198 

-.220703 

-.067590 

-.608214 

0.546450 

0.385783 

0.141075 

-.321752 

0.812732 

0.117372 

-.231162 

0.407826 

-.555128 

-.215730 

-.356564 

0.038669 

-.591894 
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TA3LE  5 

FACTOR  SCORES  OF  DATA  IN  TABLE  2 


fl 

f  2 

f  3 

1 

1.35330 

1.44370 

0.55230 

u 

0. 03500 

-1.25330 

0.74550 

3 

0.24330 

-0.34320 

1.21370 

4 

-0.23530 

0.73130 

-1.32230 

5 

-1.43420 

-0.37110 

-0.73730 

•3 

1 . 03040 

-0. 13130 

-0.33230 

7 

0.  13470 

0. 30330 

0.71440 

j 

-1.32370 

— 1 . 37030 

0.  10730 

*3 

0.20750 

-0.34500 

0.23300 

10 

- 0 . 3334  0 

-0. 17310 

0.37270 

1 1 

0.73150 

-0.24330 

1.53310 

12 

0.53350 

-0.35340 

-0. 01330 

13 

-0.22710 

-0.45330 

0 . 0333 0 

14 

—  0 . 4333  0 

-1.41430 

-0. 00730 

15 

0. 15310 

1 . 12330 

-1 . 30330 

13 

0.53430 

0.37710 

0.30350 

17 

-1. 15330 

2.23530 

0.37430 

13 

0.  03070 

0.34350 

-0.75130 

13 

0.  03330 

- 0 . 3332 0 

-0.33110 

20 

1.37330 

-  0 . 533 1 0 

- 0 . 3333 0 

21 

-0. 73450 

-1.42340 

1.15300 

22 

1.33510 

1.73750 

0. 0577 0 

0. 35340 

-0.  43170 

0. 05300 

24 

0.33530 

-1 . 04000 

1 . 233 1 0 

*jc; 

U  J 

0.23350 

1.35340 

0.  03770 

23 

0. 00430 

0 . 33  0  0 0 

0. 34070 

O'? 

• _ i 

-1.55330 

1 . 0573  0 

0.  17030 

23 

0. 03350 

-0.33730 

-0. 331 00 

23 

-1. 01370 

0.33130 

-0.2331 0 

30 

1.34030 

1.35530 

0. 23340 

31 

-0. 72540 

0.44350 

0.43450 

32 

1 . 33 120 

- 0 . 5333 0 

-0.30340 

33 

0.33450 

0.37120 

-1.23020 

34 

-1.52330 

-0. 54330 

-2.23140 

35 

-2. 02370 

0. 02750 

0.  13030 

33 

-1.24310 

-0.35300 

0.32330 

37 

-0.43100 

0.75330 

-0.51100 

33 

-0.34310 

-0.32240 

1.14530 

33 

-0.31130 

1. 13150 

-0. 13320 

4  0 

—0. 37330 

-0.53330 

1.03340 

41 

1.34010 

-1.35010 

-2.34330 

42 

0.34700 

-0. 03300 

-0.33340 

43 

0.  03150 

-0.20370 

0.31120 

44 

-0.  11050 

0.30430 

2. 13310 

45 

0.53130 

0.  33220 

0. 03220 

43 

-1.32040 

0.23010 

-1.55150 

47 

0 . 33 130 

—  1 . 32 1 3  0 

0.  0733  0 

43  . 

-0.31170 

1.33020 

-0. 00340 

43 

—  0 . 334  0  0 

-0.35130 

-1.1 0430 

50 

0.75330 

-0.30250 

1.  17300 

-  11  - 
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Since  (Aj+A^+A.^)/6  *  .78  ,  the  three  factors  explain  about  80%  of  the 

total  variability  in  the  sample.  The  communal i ties  of  the  three  factors 
are 

Variable  Communal ity 


V1 

.7593 

v5 

.9261 

v6 

.8551 

V7 

.7827 

V12 

.5184 

V14 

.8261  . 

Wo  see  that  the  three  factors  explain  more  than  75%  of  all  variables 
excluding  vJ0  .  Notice  that  the  first  factor  gives  more  weight  to  v^ 

and  v^  than  to  the  other  variables.  We  can  call  it  therefore  the  "Ship 

and  AAW  Control"  factor.  Similarly,  factor  f„  emphasizes  v,  (MW-CIC) 

and  gives  a  large  negative  weight  to  v^  .  This  factor  can  be  labeled 

"Kadar  and  Information  Communication."  The  third  factor  emphasizes  v,. 

and  dccmphasizes  v^,  v^,  v^,,  and  .  It  can  be  labeled  "Surface 

Gunnery."  In  Figure  l  we  present  the  50  simulated  vectors  represented  by 
their  first  two  factor  scores.  Such  a  representation  can  provide  a  mean¬ 
ingful  device  for  discriminating  between  "good"  and  "bad"  systems,  as 
graded  along  the  factor  scales  (f^f^.  Similar  scattergrams  can  be 

easily  provided  for  (f^,f^)  and  (f.,,f^).  As  will  be  shown  later,  such 

graphical  representation  of  the  systems  may  reveal  trends  and  clusters  of 
subsamples. 

Factor  analysis  of  multivariate  data  can  be  easily  performed  by 
employing  available  statistical  computer  packages  from  SAS,  SPSS,  or  BMDP. 
We  provide  in  Table  6  the  results  of  an  SAS  factor  analysis  procedure  per¬ 
formed  on  50  simulated  vectors  from  N(4^,4R)  ,  where  p  =  6  and  ^  is 
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TABLE  6 

FACTOR  ANALYSIS  WITH  SIX  VARIABLES, 

PRODUCED  BY  SAS  PROCEDURE  ON  50 
SIMULATED  VECTORS  FROM  N(4^,4R) 

DEFINITION  OF  VARIABLES: 

COL  1  =  ;  COL  2  =  ;  COL  3  *  ;  COL  4  *  v7  »  COL  ^  =  v12  ’  C0L  °  =  v14 

MEANS  AND  STD  DEVS 

COL  1  COL  2  COL  3  COL  4  COL  5  COL  6 
MEAN  3.91766  3.95151  4.49093  4.13207  4.32881  4.06802 

STD  DEV  1.87219  2.12153  2.02220  2.35387  2.19235  2.18691 


CORRELATION  MATRIX 


COL  l 

COL  2 

COL  3 

COL  4 

COL  5  COL  6 

COL  l 

1.00000 

0.65160 

0.13557 

0.56443 

0.28709  0.57545 

COL  2 

0.65160 

1.00000 

0.25144 

0.44143 

0.21444  0.21628 

COL  3 

0.13557 

0.25144 

1.00000 

0.12219 

0.09678  -0.08323 

COL  4 

0.56443 

0.44143 

0.12219 

1 . 00000 

0.35855  0.72271 

COL  5 

0.28709 

0.21444 

0.09678 

0.35855 

1.00000  0.23166 

COL  6 

0.57545 

0.21628 

0.08323 

0.72271 

0.23166  1.00000 

1 

2 

3 

4 

5  6 

EIGENVALUES 

2.802324 

1.156822 

0.846713 

0.685932 

0.342444  0.165765 

PORTION 

0.467 

0.193 

0.141 

0.114 

0.057  0.028 

CUM  PORTION 

0.467 

0.660 

0.801 

0.915 

0.972  1.000 

EIGENVECTORS 

1  2 


COL 

1 

0.51049 

0.04610 

COL 

2 

0.41630 

0.38367 

COL 

3 

0.12621 

0.78333 

COL 

4 

0.50975 

-0.17007 

COL 

5 

0.29683 

0.03120 

COL 

6 

0.44968 

-0.45516 

FACTOR  PATTERN 

FACTOR  1 

FACTOR  2 

COL 

1 

0.85457 

0.04958 

COL 

2 

0.69688 

0.41266 

COL 

3 

0.21127 

0.84251 

COL 

4 

0.85334 

-0.18292 

COL 

5 

0.49690 

0.03355 

COL 

6 

0.75277 

-0.48956 

FINAL  COMMUNALITY  ESTIMATES: 

COL  1  COL  2  COL  3  COL  4  COL  5  COL  6 

0.732752  0.655935  0.754461  0.761642  0.248036  0.806322 
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the  sample  correlation  matrix  of  Table  3.  The  results  obtained  are  similar 
to  those  presented  earlier.  For  the  purpose  of  simulating  the  multivariate 
data  sets  and  applying  factor  analysis  on  the  simulated  data,  we  have  found 
that  the  SAS  package  is  most  convenient.  With  the  features  available  on  SAS 
we  could  conveniently  execute  simulation  and  factor  analysis  and  other  statis¬ 
tical  procedures  in  one  program  (see  Appendix  III). 

3.2  Rotated  Factor  Analysis 

One  of  the  main  problems  of  principal  factor  analysis  is  that  the  m 
factor  score  variables  (3.1)  are  often  linear  combinations  which  ascribe  high 
relative  weight  to  many  of  the  original  variables,  and  no  immediate  (or  direct) 
interpretation  can  be  given  to  the  factor  scores.  In  order  to  obtain  factor 
scores  which  depend  on  a  small  number  of  the  original  variables  various  rota¬ 
tional  methods  have  been  developed  (see  Harman  [8]),  which  yield  vectors  of 
factors  coefficients  with  a  large  number  of  elements  close  to  zero.  We  con¬ 
sider  here  orthogonal  transformations  which  reduce  the  pattern  matrix  £  to 
a  rotated  factor  space  matrix  ^  =•  ,  where  is  an  m  x  m  orthogonal 

matrix  determined  so  that  the  column  vectors  of  ^  have  as  many  zero  entries 

as  possible.  Statistics  packages  provide  various  options  for  orthogonal  and 
oblique  rotation  of  the  factor  space.  In  Table  7  we  show  the  result  of 
EQUAMAX  of  orthogonal  rotation  of  the  factor  pattern  matrix  of  Table  6. 

This  rotation  is  designed  to  maximize  an  adjusted  fourth  moment  of  the  elements 
of  the  column  vectors  of  the  resulting  factor  pattern  matrix  ^  (see  Harman 
[8,  p.  299].  Another  commonly  applied  rotation  is  called  VARIMAX,  which  maxi¬ 
mizes  the  fourth  moment  (unadjusted)  of  each  column  of  ^  .  Both  methods  of 
rotation  are  frequently  applied  without  yielding  significant  difference  in 
the  results.  For  a  comparison  of  various  orthogonal  rotations,  including 
the  VARIMAX  and  EQUAMAX  see  Harman  [8],  Notice  that  the  factor  scores  corres¬ 
ponding  to  the  rotated  factor  analysis  are  obtained  by  multiplying  the  stan¬ 
dardized  individual  vectors  £  by  the  matrix,  .  The  "Rotated  Factor 

Pattern"  matrix  in  Table  7  is  the  matrix  ^  ^  ,  where  ^  is  the  "Factor 

Pattern"  matrix  of  Table  6  and  £  is  the  "Orthogonal  Transformation  Matrix" 
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TABLE  7 

EQUAMAX  ROTATION  OF  FACTOR  PATTERN  MATRIX 
RELATED  TO  EXAMPLE  OF  TABLE  6 


•  ROTATION  METHOD!  EQUAMAX 

ROTATED  FACTOR  PATTERN 


FACTOR 1  FACTORS 


COLl 

0. 30537 

0.29005 

COLS 

0.35114 

0.59345 

COL  3 

-0. 03650 

0.  36733 

COL  4 

0.37016 

0. 06675 

COLO 

0. 46633 

0.  17313 

0.  36074 

-0.25532 

ORTHOGONAL  TRANSFORMATION  MATRIX 

1  S 

1  0.95339  0.23377 

2  -0.23377  0.93339 


PROPORTIONAL 


CONTRIBUTIONS  TO  COMMON  VARIANCES  BY  ROTATED  FACTORS 


FACTOR 1  FACTORS 

2.369317  1.239329 


SCORING  COEFFICIENT  MATRIX 


FACTOR  1  FACTORS 


COLl 

0.23025 

0. 12763 

COLS 

0. 1 3723 

0. 41262 

COL  3 

-0. 13433 

0.71975 

COL  4 

0. 33636 

-0. 06521 

C0L5 

0. 16130 

0. 07313 

COLS 

0.37767 

-0.32357 
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of  Tabic  7.  The  "Scoring  Coefficient  Matrix"  of  Table  7  is  equal  to 


^  .  Its  two  columns  provide  the  coefficients  with  which  to  multiply 
the  individual  vectors  to  obtain  their  factor  scores.  Inspection  of 
the  scoring  coefficients  in  Table  7  shows  that  the  first  (rotated)  factor 


emphasizes  variables  v^,  v^ ,  and  v^ 


while  the  second  (rotated)  factor 


emphasizes  variables  v,  and  v. 


Variable  v 


12 


does  not  attain  a  con¬ 


siderable  weight  in  this  representation.  We  further  illustrate  the  method 
of  rotated  factor  analysis  by  performing  such  an  analysis  on  all  21  0R1 
variables.  In  Tables  8-10  we  present  the  results  of  such  an  analysis  with 
EQUAMAX  rotation.  The  data  set  consists  of  50  simulated  vectors  of  21 
components  following  the  distribution  N(5^,4^),  where  £  is  the 

correlation  matrix  of  Table  1. 


In  the  present  analysis  we  display  only  the  first  three  principal 
factors  and  their  orthogonal  rotation.  As  seen  in  Table  8,  the  first  three 
principal  factors  account  for  only  38.7%  of  the  total  variability.  In  order 
to  account  for  80%  of  the  variability  we  have  to  retain  ten  principal  fac¬ 
tors.  This  is  not  surprising,  in  light  of  the  rather  small  intercorrelations 
between  many  of  the  21  ORI  variables  (see  Table  1).  It  is  very  difficult 
to  ascribe  meaning  to  the  principal  factors  without  rotation  (see  Table  9). 
However,  after  an  EQUAMAX  rotation  we  obtain  factor-scores  coefficients 
which  can  provide  relevant  interpretation  to  the  factors  (Table  10), 
although  this  interpretation  is  different  from  the  one  obtained  by  analyzing 
six  variables  only.  Thus,  (rotated)  factor  1  emphasizes  variables 
v  ,  v]4*  vj^  and  to  some  extent  also  v?,  v^,  v^,  V2Q  •  Most  of  these 

variables  relate  to  different  aspects  of  navigation,  piloting,  and  anti¬ 
submarine  warfare.  (Rotated)  factor  2  emphasizes  v.  ,  v,,  vir>,  v.  and 

4  o  1U  11 

deemphasizes  v^q  ,  (anti-air  warfare),  engineering,  and  damage  control. 

(Rotated)  factor  3  emphasizes  v^,  v^,.,  and  v^j  ,  which  relate  to  electronic 

operations  and  setting  material  conditions.  As  explained  earlier,  not  all 
the  aspects  of  the  operational  readiness  are  represented  by  the  three 
rotated  factors.  One  needs  about  ten  rotated  factors  to  account  for  a  large 
portion  of  the  variability  in  21  variables. 
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TABLE  8 


EIGENVALUES  AND  EIGENVECTORS  OF 
CORRELATION  MATRIX  IN  TABLE  l 


EIGENVALUES 

PORTION 

CUM  PORTION 

1 

3.794984 

0.181 

0.181 

2 

2.415706 

0.115 

0.296 

3 

1.915619 

0.091 

0.387 

4 

1.732022 

0.082 

0.469 

5 

1.520571 

0.072 

0.542 

6 

1.319101 

0.063 

0.6C5 

EIGENVALUES 

PORTION 

CUM  PORTION 

7 

1.285502 

0.061 

0.666 

8 

1.212352 

0.058 

0.724 

9 

0.955083 

0.045 

0.769 

10 

0.833270 

0.040 

0.809 

11 

0.756898 

0.036 

0.845 

12 

0.633716 

0.030 

0.875 

EIGENVALUES 

PORTION 

CUM  PORTION 

13 

0.597649 

0.023 

0.903 

14 

0.460369 

0.022 

0.925 

15 

0.414631 

0.020 

0.945 

16 

0.361732 

0.017 

0.962 

17 

0.289878 

0.014 

0.976 

18 

0.193707 

0.009 

0.985 

EIGENVALUES 

PORTION 

CUM  PORTION 

19 

0.147695 

0.007 

0.992 

20 

0.115658 

0.006 

0.998 

21 

0.043856 

0.002 

1.000 

EIGENVECTORS 

1 

2 

3 

COL 

l 

0.19064 

-0.02192 

0.00061 

COL 

2 

0.31186 

-0.00371 

-0.07639 

COL 

3 

0.34164 

-0.01266 

0.00310 

COL 

4 

0.22031 

-0.23721 

0.18299 

COL 

5 

0.25561 

0.15462 

-0.00189 

COL 

6 

0.14146 

-0.26245 

0.23932 

COL 

7 

0.13558 

-0.28916 

0.06381 

COL 

8 

0.31212 

0.04438 

0.17556 

COL 

9 

0.20694 

0.22207 

0.34169 

COL 

10 

0.21507 

-0.07917 

0.27990 

COL 

11 

0.17075 

-0.25016 

0.26444 

COL 

12 

0.25605 

0.24543 

-0.01065 

COL 

13 

0.24580 

-0.17586 

-0.26415 

COL 

14 

0.22924 

-0.13077 

-0.32797 

COL 

15 

0.02812 

0.36811 

0.17200 

COL 

16 

0.31317 

0.04039 

-0.34390 

COL 

17 

0.10458 

0.29794 

-0.03943 

COL 

18 

0.23165 

0.08198 

-0.03293 

COL 

19 

0.16118 

-0.16203 

-0.16226 

COL 

20 

0.07772 

0.33043 

-0.41069 

COL 

21 

0.08802 

0.41257 

0.28378 
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TABLE  9 

COMMUNALITIES ,  FACTOR  PATTERN  MATRIX,  AND  TRANSFORMATION 
MATRIX  FOR  FACTOR  ANALYSIS  OF  TABLE  1 


SOLI 
0. 139092 
COL  3 
0.433502 
COL  15 
0.  337011 


FINAL  COMMUNAL  I  TV  ESTIMATES: 


C0L2 

COL  3 

S0L4 

C3L5  CDL6 

0.  330294 

0.443340  0. 

334266  0. 

305703  0.352044 

COL  9 

COL  1  0 

COL  1 1 

SUL  12  COL  13 

0.505296 

0.340764  0. 

395733  0. 

394545  0.437672 

COL  16 

COL  17 

COL  13 

COL  1 9  C0L2  0 

0.602679 

0.253913  0. 

221960  0. 

212452  0.609789 

FACTOR 

PATTERN 

FACTOR 1 

FACTORS 

FACTORS 

COL1 

0.37139 

-0. 03407 

0. 00035 

C0L2 

0.60752 

-0. 00577 

-0. 10572 

C0L3 

0.66553 

- 0 . 0 1 963 

0. 00429 

COL  4 

0.42917 

-0. 36369 

0.25327 

COLS 

0.49795 

0.24032 

-0.  00262 

C0L6 

0.27557 

-0. 40791 

0.33123 

C0L7 

0.26412 

-0. 44943 

0. 03332 

COL  3 

0.  6030  3 

0. 06397 

0.24299 

C0L9 

0.40313 

0.34516 

0.47292 

COL  1  0 

0.41393 

-0. 12305 

0.33740 

COL  11 

0.33264 

-0. 33332 

0 . 366  0 0 

COL  12 

0. 49331 

0.33146 

-0. 01474 

COL  13 

0. 47334 

-0.27334 

-0. 36561 

COL  14 

0.44653 

-0.20325 

-0. 45393 

COL  1 5 

0. 05473 

0.57214 

0.  23305 

COL  16 

0.61007 

0. 06273 

-0. 47597 

COL  17 

0.20372 

0. 46307 

-0. 05453 

COL  13 

0. 45127 

0. 12742 

-0. 04553 

COL  19 

0.31400 

-0.25133 

-0.22453 

COL  20 

0. 15140 

0.51357 

-0.56343 

C0L21 

0. 17147 

0.64124 

0.39277 

SOL? 
0. 279543 

COL  14 
0. 446793 
COL  21 
0.594364 


□RTHCU5DNAL  TRANSFORM AT I ON  MATRIX 
1  2 


3 


0.75536 
-0. 14324 
-0.63773 


0.43159 
-0. 61965 
0.65557 


0.49235 
0.77076 
0. 40433 


PROPORTIONAL  CONTRIBUTIONS  TO  COMMON  VARIANCES  BY  ROTATED  FACTORS 

FACTOR 1  FACT0R2  FACTORS 
3.000329  2.457702  2.663277 
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TABLE  10 

ROTATED  FACTOR  MATRIX  AND  FACTOR-SCORES  COEFFICIENTS 
FOR  FACTOR  ANALYSIS  OF  TABLE  1 


ROTATED  FACTOR  PATTERN 
FACTOR 1  FACTORS  FAC TOR 3 


COL  I 

0.S35S3 

0. 13195 

0. 15593 

COLS 

0. 53743 

0. 19547 

0. 25191 

COLS 

0. 30333 

9 . 3  0224 

0. 31425 

C0L4 

0.31753 

0.57972 

0. 02955 

COLS 

0.34343 

0. 05423 

0.42933 

COLS 

0. 05733 

0.53334 

-0. 04473 

C0L7 

0 . 3  0994 

0.45037 

-0. 13054 

COLS 

0. 39440 

0.37393 

0.  45079 

COL  9 

-0. 043 OS 

0.27014 

0.55575 

COL  10 

0. 03737 

0.51104 

0 . 253 1 1 

COL  11 

0. 07553 

0. 32443 

0 .  0 121  0 

CO- 1  s 

0. 33933 

-0. 03075 

0. 53355 

COL  13 

0 . S35S2 

0. 13 S3 3 

-0. 1 2273 

COL  14 

0.35715 

0. 02110 

-0.  12035 

COL  15 

-0. 19533 

-0. 17432 

0. 5  0421 

COL  IS 

0. 75535 

-0. 03 7 S3 

0. 15523 

COL  17 

0.  12015 

-0.23479 

0.  43514 

CO- 13 

0. 35123 

0. 03593 

0.30193 

COL  1 9 

0. 41739 

0. 14434 

-0. 13032 

COLS  0 

0.40031 

- 0 . 52553 

0.24052 

COLS1 

-0.21593 

-0. 05535 

0.  73750 

SCORING  COEFFICIENT  MATRIX 


FACTOR 1 

FACTORS 

FRCTDR3 

0. 07573 

0. 05127 

0. 03749 

0. 15555 

0.  03439 

0. 05455 

0. 13233 

0. 03220 

0. 03097 

0.  02379- 

0.23005 

-0. 00349 

0. 03530 

-0. 00591 

0. 14073 

-0. 03035 

0.24933 

-0. 02447 

0.  05073 

0. 17554 

-0. 09043 

0. 03593 

0. 13451 

0. 15219 

-0. 09333 

0. 1 1915 

2622^ 

-0. 03797 

0.21179 

0. 09533 

-0. 03173 

0.25232 

-0. 00354 

0. 03035 

• 0 • 04616 

0. 13331 

0. 23335 

.  -0.  00055 

-0. 10227 

0.25254 

-0. 05242 

-0.  10274 

-0.  10345 

-0. 05905 

0.23991 

0.27511 

-0. 10951 

-0.  00130 

0.  03033 

-0.11  429 

0 . 1 3255 

0.  0972  4 

0. 00304 

0. 03953 

0. 15275 

0. 02345 

-0. 03702 

0. 13733 

-0.30904 

0. 05351 

-0. 13595 

-0. 01 057 

0. 30975 

20  - 


T-412 


4.  Detecting  Deterioration  in 
Readiness  and  Subgroups 

We  have  seen  in  the  previous  section  that  the  readiness  of  systems 
can  be  represented  by  principal  factors  or  rotated  factors.  This  is  a  com¬ 
bined  measurement  of  readiness,  which  transforms  the  basic  ORI  scores  and 
reduces  them  to  a  small  number  of  orthogonal  factor  scores.  This  represen¬ 
tation  of  the  readiness  of  systems  is  particularly  useful  for  control  pur¬ 
poses.  Suppose  that  we  wish  to  follow  the  state  of  readiness  of  a  particular 
system.  We  can  periodically  make  observations  on  the  ORI  variables  and  pre¬ 
sent  the  corresponding  factor  scores  on  the  scattergrams  similar  to  the  one 
in  Figure  1.  Significant  deterioration  in  readiness  will  be  detected  by 
the  location  of  these  points  in  the  scattergram.  Moreover,  if  a  whole  group 
of  points  cluster  on  the  scattergram  on  the  negative  side  of  a  factor  there 
may  be  an  indication  that  this  group  originates  from  a  different  population 
and  further  analysis  should  follow.  Such  a  case  is  illustrated  in  Figure  2, 
in  which  the  factor  scores  obtained  by  a  rotated  factor  analysis  of  the  6 
ORI  variables,  when  the  sample  of  50  systems  consisted  of  25  units  from  the 
distribution  N(5^,4£)  and  25  units  from  the  distribution  NQ,4£).  The 

points  in  Figure  2  corresponding  to  the  units  in  the  first  subsample  are 
labeled  "1"  and  the  others  are  labeled  "2".  It  is  seen  that  most  of  the 
second  subsample  points  are  concentrated  at  the  negative  part  of  f^  . 

There  is  a  strong  indication  of  a  significant  difference  between  the  two 
subsamples.  The  capability  of  rotated  factor  analysis  to  separate  such  sub¬ 
samples  in  the  new  factor  space  is  not  surprising.  It  can  be  given  precise 
algebraic  and  geometric  interpretations.  We  do  not  dwell  on  this  here  to 
any  further  degree  but  only  remark  that  if  such  separation  of  two  natural 
subsamples  (as  two  different  subfleets)  is  indicated  then  one  should  rein¬ 
force  the  analysis  by  performing  another  method  of  multivariate  analysis, 
which  is  designed  for  discrimination  between  subgroups  and  classification 
of  the  individual  units  to  various  readiness  groups  according  to  their  dis¬ 
tances  from  the  centroids  of  these  groups.  For  the  theory  and  explanation 
of  these  methods  see  Tatsuoka  [12],  Afifi  and  Azen  [1],  and  Van  De  Geer  [13]. 
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APPENDIX  I 

SIMULATING  MULTINORMAL  VECTORS 

In  the  present  appendix  we  discuss  use  of  the  SAS  computer  package 
for  simulation  of  p-dimensional  multinormal  vectors.  PROCEDURE  MATRIX  of 
SAS  provides  the  possibility  of  computing  the  eigenvalues  and  an  ortho¬ 
normal  matrix  of  eigenvectors  of  the  correlation  matrix  ^  .  This  means 
that  the  simulation  can  be  based  on  the  following  result  from  the  theory 
of  multinormal  distributions. 

Let  be  an  orthonormal  matrix  of  the  eigenvectors  of  a  correlation 
matrix,  ,  and  let  |  be  a  diagonal  matrix  whose  elements  are  the  eigen¬ 
values  of  £  .  Assume  that  is  of  full  rank.  We  can  write  ^  T  . 

1/2 

Accordingly,  let  and  let  ^  be  a  p-dimensional  vector  of 

independent  standard  normal  random  variables,  that  is,  £  %  >  then 

the  distribution  of  %  =  is  N(£,£)  .  Finally,  ^  +  2%  is  dis¬ 

tributed  like  N(u^,4j^)  .  The  SAS  program  given  in  Table  1,  Appendix  III 
simulates  50  independent  21-dimensional  multinormal  vectors  with  mean  vector 
5,|  and  covariance  matrix  4,g  ,  where  ^  is  the  correlation  matrix  of  Table  1 
and  it  then  performs  the  rotated  factor  analysis.  The  program  is  based  on 
the  following  algorithm: 

Step  1.  Read  jj.  ; 

Step  2.  Generate  N  *  50  vectors  of  21  independent 
standard  normal  random  variables; 

Step  3.  Determine  the  eigenvalues  and  normalized 
eignvectors  of  ; 

Step  4.  Arrange  the  eigenvalues  in  a  diagonal  matrix  J) 
and  the  eigenvectors  in  a  matrix  ; 

1/2 

Step  5.  Determine  £  =  ] £  and  £  *  ; 

Step  6.  Arrange  the  data  generated  in  Step  2  in  a  50  x  21 
matrix  %  ; 
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Step  7.  Make  the  transformation  W  =  2*Y»C  : 

Step  8.  Determine  the  matrix  jjl  =  5*^  ,  where  ^  is 
a  50  x  21  matrix  of  L's; 

Step  9.  Compute  • 

The  matrix  consists  of  50  i.i.d.  row  vectors,  each  of  which  is  distributed 
like  N(5^4£)  . 

If  a  SAS  package  is  not  available  one  can  perform  the  simulation  by 
another  method  which  does  not  require  the  determination  of  eigenvalues  and 
eigenvectors  but  only  the  solution  of  linear  equations.  A  FORTRAN  program  of 
such  a  procedure,  based  on  a  recursive  algorithm,  is  given  in  our  report  [4]. 
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PRINCIPAL  AND  ROTATED  FACTOR  ANALYSIS 


be  a  standard  multivariate  normal  vector.  The 


distribution  of 


be  the  eigenvalues  of  ^ 
with  length  ||£||  =  1  , 


.  We  wish  to  determine  a  vector  (functional)  £ 
which  maximizes  the  variance  of  £Tu  .  The 


ferentiating  f(^»A)  with  respect  to  ^  yields  the  eigenstructure  equation 


Notice  that 


Thus,  £ 


is  an  eigenvector  of  ,  of  unit  length,  corresponding 
value  of  R  ,  namely  to  A.  which  is  the  variance  of 


Similarly,  let  ' 
corresponding  to  A 


be  the  eigenvectors  of  unit  length 
Notice  that  the  variance  of 


Indeed,  if 


decomposition  of  ^  is 
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Furthermore,  for  any  i  j*  j 


■ },  ^<i)T«<J)  ■ 0 

k=l 


(A. 2.4) 


Let  ^  be  an  orthogonal  matrix  with  columns  which 

are  the  eigenvectors  of  £  .  The  distribution  of 


l  =  (A-2-5) 

is  like  that  of  ;  where  ^  =  diag(Ar . . .  ,A  )  .  Indeed, 

T  T  T 

^  J^jJ)  .  But  J^S  ^  =  ^  *  The  comPonents  of  £  are  called  the 

principal  factor  scores,  corresponding  to  ij  .  The  orthogonal  transformation 

°f  A!  »  8iven  by  (A. 2. 5)  yields  independent  standard  normal  random  variables. 

P 

Since  trace  £  =  trace  A  =  E  A  =  p  ,  the  ratio  A./p  (i=l,...,p) 

i=l  1  1 

is  the  proportion  of  the  total  variance  of  accounted  for  by 

f^  (i=l,...,p)  .  If  we  choose  only  the  first  m  (l<m<p)  eigenvectors  of 

£  ,  corresponding  to  \1  >...>  Am  ,  and  define  £  =  [£  (1) , . . . ,£(m) ]  , 

A(m)  =  diag(^r---»^m)  .  the  transformation  %  yields  the  first 

m  components  of  ^  .  The  concepts  of  communality  and  the  nature  of  rotated 
factor  analysis  was  explained  in  Section  3. 


* 
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COMPUTER  PROGRAMS 

In  the  present  appendix  we  present  two  SAS  programs.  The  program 
in  Table  III.l  performs  simulation  of  50  independent  multinormal  vectors 
of  21  components,  having  the  common  distribution  N(5^,4J^)  .  A  rotated 
factor  analysis  is  then  performed  on  the  simulated  data  set  with  an  EQUAMAX 
rotation.  In  Table  III. 2  we  present  a  program  according  to  which  a  scatter- 
gram  of  the  factor  scores,  corresponding  to  the  simulated  data  set,  can  be 
obtained.  This  program  is  designed  to  present  the  scattergram  of  f2 
versus  f^  in  the  case  of  two  subsamples  of  size  n  =  25  from  N(4^,4j^) 
N(2.5^.,4J^)  ,  respectibely .  Figure  2  was  obtained  by  a  similar  program,  with 
two  subsamples  from  N(5^.,4^)  and  N(^,4j^)  .  The  points  from  the  two 
subsamples  are  labeled  "1"  and  "2",  respectively.  The  first  part  of  the 
program  simulates  50  6-dimensional  normal  vectors  from  these  two  distributions. 
The  factor  weights,  to  obtain  f^  and  f ^  ate  read  into  the  program  as  a 
data  set  ("DATA  FSCORE") .  These  weights  can  be  obtained  by  performing  the 
factor  analysis  program  in  Table  III.l. 
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TABLE  1 1 1.1 

SAS  PROGRAM  FOR  SIMULATING  'jU  N(5^,4£)  VECTORS  AM) 
PERFORMING  ROTATED  FACTOR  ANALYSIS  ON  21  VARIABLES 

□ptidns  ls*3o; 

DATA  cqrri 
.INPUT  Cl -call 

TITLE  FACTOR  ANALYSIS  FOR  READINESS? 

CARDS! 


i. 

00 

0. 

17 

0. 

27 

0. 

04 

0. 

33 

0. 

31 

0. 

45 

0. 

10 

0. 

06 

0. 

03 

0. 

03 

0. 

32 

0. 

36 

0. 

36 

0. 

17 

0. 

16 

0. 

07 

3. 

12 

0. 

03 

0. 

01 

0. 

10 

0. 

17 

1. 

00 

0. 

22 

0. 

06 

0. 

13 

0. 

03 

0. 

14 

0. 

22 

0. 

01 

0. 

20 

0. 

04 

0. 

12 

0. 

14 

0. 

26 

0. 

01 

0. 

07 

0. 

02 

0. 

09 

0. 

09 

0. 

03 

0. 

13 

0. 

97 

k_  1 

0. 

22 

i. 

00 

0. 

06 

0. 

13 

0. 

03 

0. 

14 

0. 

22 

0. 

01 

0. 

20 

0. 

04 

0. 

12 

0. 

14 

0. 

26 

0. 

01 

0. 

07 

0. 

02 

0. 

09 

0. 

09 

0. 

03 

0. 

13 

0. 

04 

0. 

06 

0. 

06 

1. 

00 

0. 

23 

0. 

23 

0. 

13 

0. 

11 

0. 

13 

0. 

10 

0. 

03 

0. 

12 

0. 

35 

0. 

03 

0. 

19 

0. 

25 

0. 

03 

0. 

10 

0. 

03 

0. 

04 

0. 

06 

0. 

33 

0. 

13 

0. 

13 

0. 

23 

1. 

00 

0. 

34 

0. 

05 

0. 

03 

0. 

04 

0. 

01 

0. 

12 

0. 

07 

0. 

03 

3. 

04 

0. 

04 

0. 

25 

0. 

10 

0. 

20 

0. 

03 

0. 

13 

0. 

02 

0. 

31 

0. 

03 

0. 

03 

0. 

23 

0. 

34 

1. 

00 

0. 

29 

0. 

13 

0. 

13 

0. 

03 

0. 

05 

0. 

03 

0. 

22 

0. 

15 

0. 

06 

0. 

17 

0. 

12 

0. 

13 

0. 

04 

0. 

03 

0. 

10 

0. 

45 

0. 

14 

0. 

14 

0. 

13 

0. 

05 

0. 

29 

1. 

00 

0. 

17 

0. 

17 

0. 

03 

0. 

09 

0. 

14 

0. 

17 

0. 

11 

0. 

17 

0. 

14 

0. 

01 

0. 

31 

0. 

01 

0. 

12 

0. 

03 

3. 

10 

0. 

22 

0. 

22 

0. 

11 

0. 

03 

0. 

13 

0. 

17 

1. 

00 

0. 

09 

0. 

1  0 

0. 

00 

0. 

21 

0. 

4  0 

'3. 

01 

0. 

03 

0. 

10 

0. 

17 

0. 

23 

0. 

12 

0. 

07 

0. 

03 

3. 

06 

0. 

01 

0. 

01 

0. 

13 

0. 

04 

0. 

13 

0. 

17 

0. 

09 

1. 

00 

0. 

13 

0. 
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TABLE  III.l  (Cont'd) 


DRTR  RRND? 

XEEP  X1-X21? 

loopjn+i; 

X 1  =NORMRL  < 0  -  ;  X2=N0RMRL  <  0 >  5  X3-N0RMRL  <  0>  5  X4=N0RMRL  <  CP  5  X5=N0RMRL  (  0)  5  X6=*N0RMRL  <  0  >  5 
X?=N0RMRL  <  0>  ?  X3=N0RMRL  <  0>  5  X9=N0RMRL  <  Q>  5  X 1  Q-NDRMRL  <  O  '  5 

XI 1 -NORMAL  <0>  5  X12=NORMRL  <0  :-  5  X13=NORMAL  <0>  ?  X 1 4=N0RMRL  ( 0>  ?  X15=NQRMRL  <0>  5 
X 1  o-NORMAL  <  0>  i  X 1 7=N0RMAL  (  0  :•  ?  X 1 3=N0RMRL  <  CP  5  X 1 3=N0RMAL  <  0  >  ?  X2  0=NORMAL  <  0>  ? 

X2 1 -NORMAL  <  0  > 5 
□UTPIJT? 

IF  N  <  50  THEN  GO  TO  LOOP 5 
PROC  MATRIX? 

FETCH  R  HRTA=CORR? 

EIGEN  M  E  R5 
d=diag<mj ? 

S»SQRT <D>  ; 

3=E»3? 

FETCH  Y  DATAGRAMS? 

♦b  ■- ; 

w=2*w? 

I A- 1 ! 5  0  ?  W 1 =W  < I A  j ♦>  ?  MW 1 =  J  <5  0  >  2 1 j  5>  5 

w=wi »nui ; 

OUTPUT  W  OUT =D  I  SC  <KEEP=COL  1 -C0L2 1  >  5 

PROC  F RC TOR  NFACT=3  OUT=FACTl  METHOD=PRlN  RDTATE=EQUAMAX  EIGENVECTORS  SCORE? 
7RR  C0L1-C0L21? 
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TABLE  II I. 2 

SI  Mil  I  .AT  INC  MULT  I  NORMAL  VECTORS  WJT11  SIX  COMPONENTS  FROM  TWO 
DISTRIBUTIONS  AND  PLOTTING  THEIR  FACTOR  SCORES 


3PTI0NS 
DATA  :3RP? 

INPUT  C1-C6*# 

TITLE  FACTOR  ANALYSIS  FOR  READINESS? 

cards; 

1.0000  0.5354  0.1339  0.5532  0.3390  0.3195 
0 . 535 4  1.0000  0.  1 532  0 .  3623  0 . 2334  0.1231 
0.  1339  0.  1532  1.0000  0.2374  0.  1927  -.  0234 
0 . 5532  0 . 3623  0 . 2374  1.0000  0 . 492 1  0.6177 
0 .3390  0.2334  0.  1 927  0 . 492 1  1.  000  0  0 . 2399 
0.5195  0.1231  -.  0234  0.6177  0.2399  1.0000 
DATA  RAND5 
KEEP  X1-X65 
_33P:N*i; 


:<  i  -normal  <  o  • ;  X2 -normal.  <  o  • ;  x  3  -normal  o  • ;  x  4  -normal  <  o  < ;  X5  -normal  <  o> ;  ms-normal  <  o » 


□■jtput; 

IF  N  <  50  THEN  S3  TO  L33P; 

DATA  f score; INPUT  F1-F2;CARDS5 

0.49979  -.014726 

0.52570  -.024939 

-.21200  0.43125 

0.14143  .27133 

-.  13526  0.50603 

0.17949  0.17796 

proc  matrix; 

FETCH  R  DATA=CCJRR; 

EI5EN  M  E  R; 

D-DIAG<M>  ; 

S=SQRT  <D> ? 

b=e>s; 

FETCH  v  data=rand; 

W=Y*By; 

W=2 ♦W ; 

I A-l :  25?  W  1  =1.1  CIA*  ♦>  ;  MW  1  =  J  <25 >  6 >  4>  ? 
wi=wi>mwi; 

I  B=26 : 5  0 ;  W2=W  <1 B  t  ♦>  ;  MW2- J  <25 . 6  >  2 . 5>  *>  W2-W2+MW2 ; 
w=wi/'W2; 

i  D= i  <5  o:> ;  W3=j  <5  o  >  5  o » i  > ;  w  3 = .  02H13 ; 

b=id-w3;d=b*w;s=w"»d;v=diag(S>  ;s=sqrt<7>  ;ij=inv(S'  ;z=d*u; 

FETCH  F  DATA=FSC3RE;SC=Z»F? 

OUTPUT  SC  3UT-SC0R (KEEP-COL 1 -C0L2> ? 

DATA  SC3R;SET  SCOR? 

IF  _N_  LE  25  THEN  SP=1 ? IF  _N_  GT  25  THEN  GP-2? 

PROC  PLOT? PLOT  C3L2+C3L 1 -GP ? 
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